Azimuthal distributions of radial momentum and velocity in relativistic heavy ion 

collisions 
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Azimuthal distributions of radial (transverse) momentum, mean radial momentum, and mean 
radial velocity of final state particles are suggested for relativistic heavy ion collisions. Using trans- 
port model AMPT with string melting, these distributions for Au + Au collisions at 200 GeV are 
presented and studied. It is demonstrated that the distribution of total radial momentum is more 
sensitive to the anisotropic expansion, as the anisotropics of final state particles and their associated 
transverse momentums are both counted in the measure. The mean radial velocity distribution is 
compared with the radial flow velocity. The thermal motion contributes an isotropic constant to 
mean radial velocity. 

PACS numbers: 



I. INTRODUCTION 

One of the main goals of current relativistic heavy ion 
collisions is to understand the properties of a new formed 
matter — the quark-gluon plasma (QGP) It is well 
known that one important character of the formed mat- 
ter is anisotropic collective flow. In non-central collisions, 
the overlap area of two incident nuclei is an almond shape 
in the transverse coordinate plane Q. This initial geo- 
metric asymmetry leads larger density gradient along the 
short axis. It in turn pushes the formed system to expand 
anisotropically, i.e., large collective flow velocity in short 
side direction, which is perpendicular to the anisotropy in 
coordinate space. Therefore, the measure of anisotropic 
distribution of final state particles shouldprovide valu- 
able information of the system evolution [J Q . 

Conventionally, the azimuthal distribution of multi- 
plicity of final state particles is presented. Its anisotropy 
is quantified by the coefficients of Fourier expansion of 
the distribution Q 
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where <j> is the azimuthal angle between the transverse 
momentum of the particle and the reaction plane. The 
Fourier coefficients is evaluated by, 



v n (N) = (cos(n^)), 



(2) 



where (...) is an average over all particles in all events, 
and v n (N) refers to the anisotropy coefficient of az- 
imuthal multiplicity distribution. The second harmonic 
coefficient V2(N) is the so-called elliptic flow parameter. 
It presents the anisotropy of the colliding system and 
has the biggest ellipticity at high energy heavy ion col- 
lisions [H, Besides, the azimuthal asymmetry distribu- 
tion of energy loss and its coefficients of Fourier expan- 
sion are also studied 01 . 



However, the multiplicity distribution only counts the 
number of particles emission in a certain azimuthal an- 
gle. The expansion of the system results in not only the 
anisotropy of multiplicity distribution but also their as- 
sociate radial (transverse) momentum. The total radial 
momenta at a given azimuthal angle is the combination 
of them. Therefore, the azimuthal distribution of radial 
momentum should be a more sensitive measure of the 
anisotropic expansion, which has not been directly ex- 
plored before. 

In addition to the radial momentum, the radial flow 
velocity is another interesting and important quantity. 
It directly relates to the equation of state Q and shear 
viscous interactions. For an ideal flow, the radial flow 
velocity is isotropy. While, if there are shear interactions, 
the radial flow velocity will be different from layer to 
layer. In hydrodynamics, the shear viscous interactions 
are supposed to be proportional to the gradient of flow 
velocity @. The proportional constant is defined as shear 
viscosity. The gradient of radial flow velocity along the 
azimuthal direction is directly related to the shear viscous 
interactions. 

Theoretically, the radial flow velocity is a parameter in 
model calculations. It is usually obtained by fitting the 
spectrum of transverse momentum [lOj . Recently, it is 
further suggested to extract the radial flow velocity from 
photon and dilepton spectrum fill ]. 

Experimentally, only the radial velocity of final state 
particles (v) is measurable. It should be a combination of 
velocities of flow ("Uflow) and the random thermal motion 
(^th) [13 |. How to abstract the random thermal motion 
from the radial velocity of final state particles and get the 
radial flow velocity are not clear. This is why the radial 
velocity of final state particles has not been explored in 
a long period. It is interesting to see how the radial 
velocity of final state particle relates to the radial flow 
velocity. Therefore, we further suggest the measurement 
of the azimuthal distribution of mean radial velocity of 
final state particles. 
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In the second session of the paper, we will give the 
definitions of suggested azimuthal distributions of radial 
momentum and velocity of final state particles, and the 
corresponding anisotropic parameters. In the third ses- 
sion, using the samples generated by AMPT with string 
melting model, we show the azimuthal distributions of 
radial momentum and the centrality dependence of its 
anisotropic parameters. The results are compared with 
those of the corresponding azimuthal multiplicity distri- 
bution. In the fourth session, the azimuthal distributions 
of mean radial velocity at different ccntralitics arc pre- 
sented, and compared with those given by the anisotropic 
blast- wave model fl3l - fl5j . Finally, the summary and con- 
clusions are presented. 



II. AZIMUTHAL DISTRIBUTIONS OF RADIAL 
MOMENTUM AND VELOCITY 

As indicated, the initial anisotropy in coordinate space 
in non-central collisions makes the formed system expand 
in a perpendicular almond shape in momentum space. 
The final state particles move outward anisotropically. 
Both the particle density and the associated momentum 
behaves anisotropically during the expansion. The dis- 
tribution of total transverse momentum at the different 
azimuthal directions should be a good measurement for 
both of these two effects. The total transverse momen- 
tum in the mth azimuthal bin can be defined as 
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where pt,; is the transverse momentum of the ith particle, 
N m is the total number of particles, and (. . .) denotes the 
average over all events. 

In order to see the contributions of radial momentum 
in particular, the mean radial momentum in the mth 
azimuthal bin can be defined accordingly as, 



«Pt(0 m )» 



N, 



1 Movent / 1 N m \ 

— E UrE^(^) ■ (4) 



Here, the averages ((. . .)) arc over all particles in mth 
angle bin and all events. It records only the contribu- 
tions from the transverse momentum of final particles, 
the multiplicity effect is canceled by the average over all 
particles. 

The anisotropic parameters of all those azimuthal dis- 
tributions can be directly obtained from their Fourier 
expansions, respectively, 

^oc l + f>„«P t »cos(r^), (5) 



n=l 



and 



^g^-oc l + f>„(«p t »)eos(n«£). (6) 



and ^Bp- are the azimuthal distribution functions 
of total radial momentum and mean radial momentum. 
v n ((Pt)) a- n d v n (((pt)}) are their anisotropic parameters, 
respectively. 

Considering the rclativistic effect, the transverse (ra- 
dial) velocity of the ith particle can be written as, 



Pt,i 



Pi. 



(7) 



Pi 



where pt,i and $ are the transverse momentum and 
mass of the ith particle, respectively. mo,i is the mass of 
ith particle in rest frame. The radial velocity fluctuates 
from particle to particle. In a given azimuthal direction, 
the mean radial velocity can be considered as a good 
approximation. Analogously, the azimuthal distribution 
of mean radial velocity can be defined as 



N, 



— E (^E U M(0m)J 



(8) 



Here, the average is over all the particles in mth bin and 
events. 

The behavior of those suggested observables should 
provide more information in anisotropic expansion. In 
the following, as a demonstration, we use the generated 
sample of AMPT with string melting [Tj| [Tjj • A partonic 
phase is implemented in the model and the elliptic flow 
data from RHIC are well reproduced by the model (l8j |. 
For Au+Au at ^s^n = 200 GeV, about 1.6 millions 



(3) 

minimum bias events are generated. 



n=l 



III. AZIMUTHAL DISTRIBUTIONS OF 
RADIAL MOMENTUM IN AMPT MODEL 

The azimuthal distributions of radial momentum, 
mean radial momentum, and multiplicity are presented 
in Fig. 1(a), (b) and (c), respectively. Error is statistical 
only and smaller than the size of the points. The parti- 
cles within rapidity range y £ [—5, 5] are counted. These 
cases are kept in all the following figures. 

We can see from Fig. 1 that all the observables as a 
function of azimuthal angle show the anisotropic shape, 
cos(20). It is the same as multiplicity distribution, the 
biggest anisotropy of mean radial momentum distribu- 
tion appears in in-plan direction as shown in Fig. 1(b). 
It indicates that not only the particle density, but also 
the associated pt are larger in in-plane direction. It is 
interesting to see if the data at RHIC show the same 
character as the model. 

In order to compare the anisotropy effects of these 
three distributions qualitatively, the centrality depen- 
dence of the corresponding anisotropic parameter, 1)2, is 
presented in Fig. 2. The anisotropy parameter V2 from 
different measurements shows the similar centrality de- 
pendencies. At each centralities, the anisotropy parame- 
ter of multiplicity distribution, v-i (N) , is larger than that 
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FIG. 1: The azimuthal distributions of radial momentum (a), 
mean radial momentum (b), and multiplicity (c) for the sam- 
ple of Au+Au collisions at ^snn = 200 GeV generated by 
AMPT with string melting. 
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FIG. 2: (Color online)The centrality dependence of elliptic 
flow parameters deduced from azimuthal distributions of ra- 
dial momentum (solid red stars), mean radial momentum 
(solid black cycles), and multiplicity (solid blue triangles) for 
the sample of Au+Au collisions at ^/snn = 200 GeV gener- 
ated by AMPT with string melting. 



FIG. 3: (Color online) p t dependence of anisotropic parame- 
ter of azimuthal distributions of radial momentum (red solid 
stars), and multiplicity (black triangles) for the sample of 
Au+Au collisions at ^/sWn = 200 GeV generated by AMPT 
with string melting. 
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FIG. 4: (Color online) The azimuthal distributions of mean 
radial velocities of minimum bias sample (a) and the samples 
of three different centralities (b), and the azimuthal gradients 
of (b) in (c). 



of the mean radial momentum distribution ^(((p*))). 
The anisotropy parameter of radial momentum, V2{(Pt}), 
is the largest one among the three variables. It confirms 
that the anisotropy of radial momentum distribution in- 
cluding the contributions from a number of particles and 
their associated transverse momentum. Therefore, the 
azimuthal distribution of radial momentum gives a full 
count of anisotropic expansion. 

As we know, the anisotropy parameters vi also depend 
on p t , and it increases with p t when p t < 2 GeV/c fl9| . 
The pt dependence of the anisotropy parameters of radial 
momentum and multiplicity distributions are presented 
in Fig. 3. The anisotropy parameter increases with p t 
when p t < 2 GeV/c, the same as the data shown. We 
can also see that the vi slightly decreases with p t when 
Pt > 2 GeV 7c, and it may be contributed by the hard 
components [20[ • At a fixed pt bin, the anisotropy of the 
radial momentum is almost the same as that of multiplic- 
ity. This is because the pt of all particles in a small given 
Pt bin are almost the same. The anisotropy of radial 



momentum is dominated by that of multiplicity. 

IV. AZIMUTHAL DISTRIBUTIONS OF 
RADIAL VELOCITY IN AMPT MODEL 

The azimuthal distribution of mean radial velocity is 
presented in Fig. 4(a). It is a period function and can be 
well fitted by 

{( v t)) = Vo + V a cos(2<£). (9) 

It is the same mode as the flow velocity, 

/3 = /3 o + /? a cos(20), (10) 

which is usually assumed in blast-wave model in counting 
the anisotropic expansion (l~3l . Il5j . 

In order to see the contribution of the random thermal 
motion, the mean radial velocity of final state particles 
in three typical centralities are presented in Fig. 4(b). 
In mid-central (30% - 40%) and peripheral (60% - 70%) 
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collisions, the mean radial velocities are anisotropy, while 
it becomes approximately an isotropy constant in central 
(0—5%) collisions. This suggests that the interactions be- 
tween azimuthal layers are negligible in central collisions, 
which is consistent with the expectations of viscous hy- 
drodynamics [2l|, [22[ . It also shows that the thermal mo- 
tions only contribute an isotropic constant to the mean 
radial velocity. 

As we know, for a system with a fixed temperature, 
the lighter particle has higher thermal velocity. In order 
to test if the Vq is mainly caused by thermal motion, the 
mean radial velocities of three different particles and their 
corresponding fitting parameters are presented in Fig. 5. 
Indeed, the lightest pion has the highest Vq, while the 
heaviest proton has the lowest one. 

To see the anisotropy effect alone, we can calculate the 
gradient of mean radial velocity along the azimuthal di- 
rection. In the case, the constant part of the mean radial 
velocity is canceled. Fig. 4(c) shows the corresponding 
gradients of Fig. 4(b). In central collisions, it is approxi- 
mately zero. The amplitudes in mid-central collisions are 
larger than those in peripheral collisions. These results 
show that there is almost no gradient of mean radial ve- 
locity in central collisions and it becomes the largest in 
mid-central collisions. 

Conventionally, the parameters of flow velocity 
Eq. (|10[) . (3q and f3 a , are obtained by fitting the spectra of 
produced particles. Here, we choose the spectra of pion, 
proton and kaon from AMPT string melting and get, 
[3 = 0.35 + O.O4cos(20). Due to thermal motion, the fto 
is not directly comparable with Vq. However, V a ~ 0.01 
from corresponding mean radial velocity may be a good 
approximation of flow velocity estimated by blast-wave 
model, where (3 a ~ 0.04. 

Certainly, the flow velocities obtained from directly 
measured radial velocity and from the spectrum fitting 
based on blast-wave model should be better compared 
by experimental data sample, where the spectrum is pre- 
cisely presented. The comparison of these two methods 
will lead to a better understanding of the flow velocity. 



V. SUMMARY AND CONCLUSIONS 

In the paper, we suggest the studies for azimuthal dis- 
tributions of radial momentum, mean radial momentum, 
and mean radial velocity in relativistic heavy heavy ion 
collisions. 

Using the sample of Au + Au collisions at ^snn = 
200 GeV produced by a multiphase transport model 
(AMPT) , we find that the azimuthal distribution of ra- 



dial transverse momentum indeed counts the anisotropy 
of final state particles and their associated transverse mo- 
menta. Thus it presents a full description of anisotropic 
expansion at various centralities. Only in small pt bin, 
the azimuthal distribution of radial momentum shows the 
same anisotropy as that of the multiplicity distribution. 

The azimuthal distribution of mean radial velocity is 
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FIG. 5: (Color online)The azimuthal distributions of mean ra- 
dial velocity of charged pion, charged kaon and (anti)proton. 
The lines are fitted by Eq. (8), and the corresponding fitted 
parameters are listed. The errors of the parameters are less 
than 1% relative values. 



shown to be the same mode as flow velocity which is usu- 
ally assumed in generalized blast-wave model. Its cen- 
trality dependency indicates that thermal motion only 
contributes an isotropic constant to mean radial velocity. 
Its particle mass dependency further shows that the mass 
ordering of isotropic mean radial velocity is the same as 
thermal motion. The anisotropic mean radial velocity 
is approximated to flow velocity, which is obtained from 
fitting the spectrum of corresponding particles based on 
Blastwave model. 

Therefore, it is interesting to measure the azimuthal 
distributions of radial momentum, mean radial momen- 
tum, and mean radial velocity in current relativistic 
heavy ion collisions. 
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